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We present calculations of ground state properties of spherical, doubly closed-shell nuclei from 
16 O to 208 Pb employing the techniques of many-body perturbation theory using a separable density 
dependent monopole interaction. The model gives results in Hartree-Fock order which are of similar 
quality to other effective density-dependent interactions. In addition, second and third order per- 
turbation corrections to the binding energy are calculated and are found to contribute small, but 
non-negligible corrections beyond the mean-field result. The perturbation series converges quickly, 
suggesting that this method may be used to calculate fully correlated wavefunctions with only sec- 
(2JT)| ond or third order perturbation theory. We discuss the quality of the results and suggest possible 

^ ■ methods of improvement. 
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PACS numbers: 21.10.Dr, 21.10.Ft, 21.30. Fe, 21.60.Ev, 21.60.Jz 



I. INTRODUCTION 
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The central problem of nuclear structure theory is the solution of the many-body Schrodinger equation (MBSE). 
For Hamiltionans of interest in the nuclear case, an analytic solution is impossible, and one is compelled to use some 
approximation, either in the numerical solution of the equation or the specification of the Hamiltonian, or both. 

Approaching the problem with the aim of using as realistic a representation of the potential as possible usually means 
fitting a combination of a meson exchange and phenomenological interaction to low-energy nucleon-nucleon scattering 
data and properties of few-body systems. To get good agreement with experiment both two- and three-body forces 
seem to be necessary. Recent examples of such potentials include the Bonn [Q , the Argonne two-body || with Urbanna 
3-body Nijmegen j| and Moscow § potentials, the last of which also incorporates quark degrees of freedom. These 
forces share the property of having a hard repulsive core which is a natural consequence of meson-exchange. It is this 
hard core which presents the difficulty in solving the MBSE. For instance, Hartree-Fock (HF) mean-field calculations 
with such interactions result in unbound nuclei. Treating corrections beyond the HF approximation order-by-order 
in perturbation theory is also unsuccessful since the interactions used are non-perturbative. One has to solve the 
full MBSE numerically in as exact a way as possible using techniques such as Variational Monte-Carlo | B, Green's 
Function Monte Carlo (]]] , the coupled-cluster method , and the Fermion Hyper-netted chain model po[ . Using 
effective interactions derived from realistic potentials, no-core shell-model calculations have been made in light nuclei 
p[ and heavier nuclei close to closed shells have been treated p3| . 

The computational difficulty of performing numerically exact solutions of the MBSE has limited the techniques 
to light nuclei, for instance A = 8 results have been published recently using the Argonne ul8 and Urbanna IX 
potentials in the GFMC framework (l2) . In this work, it is seen that although the lightest nuclei are reproduced very 
well, the quantitative comparison of theory to data gets worse as A increases. This may be due to the necessarily 
phenomenological nature of the three-body potential, a problem which may be overcome with re-fitting. On the other 
hand, it is not obvious that higher-body forces will not prove necessary or that the concept of a bare interaction 
between nucleons is valid for small distances. 

Attempts were made in the late sixties primarily by the Kerman group at MIT to parameterize the NN interaction 
in such a way that it is weak in the sense of being perturbative. Such a weak interaction allows one to perform Hartree- 
Fock calculations to obtain a reasonable approximation to the full wavefunction and then to calculate corrections in 
perturbation theory. While this technique seems very attractive, the result obtained were only moderately successful 
at reproducing experimental data |p^-|l8|] , a fact which was presumed to be due to inadequacies in the potentials used. 
The efficacy of developing a suitable interaction when similar, though more complicated, techniques were available for 
realistic interactions has been questioned [ fl9| and no better interaction was developed. Separable parameterizations, 
particularly the quadrupole-quadrupole interaction [GQEll , have retained currency, as residual interactions E) . Even 
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when the interactions are too strong for regular perturbation theory, separable interactions requiring solution of 
Bruckner Hartree-Fock equations have proved fruitful [^2| because of their simplicity. 

On the other hand, interactions have been developed which are not intended for use in the full MBSE, but rather 
to give good results with a Hartree-Fock calculation alone. Good quantitative success came with the zero-range 
density-dependent force of Ehlers and Moszkowski p3] and Skyrme's interaction |23], used in HF calculations by 
Vautherin and Brink [^5| and subsequently by many others, and also Gogny's finite-range interaction [p6[ . Skyrme's 
interaction has been particularly successful, in part due to its simple form, that of a delta function, which leads to 
easy calculation, even of the exchange part of the force. This computational simplicity has allowed extensive study 
of the prope rties of nuclei to be made with the Skyrme interaction across the entire range of nuclei in the periodic 
tab le Ip7| , ff8| , [50t . Related somewhat to the Skyrme-Hartree-Fock model is the Relativistic Mean Field (RMF) approach 
p8| , p9| , which also gives single-particle motion in a mean field, but as a solution to the Dirac equation as opposed to 
the Schrodinger equation. The RMF approach has some nice features such as the natural occurrence of the spin-orbit 
splitting without recourse to an assumed spin-orbit interaction. 

These mean-field models are inherently single-particle in nature. The forces used are not intended to be used in the 
MBSE, nor are explicit corrections beyond the mean-field part of the framework, although Skyrme's interaction can be 
considered as a phenomenological G-matrix equivalent |30[ , and in that sense includes a subset of possible correlation 
effects in the mean-field. Although appropriate for use in mean-field calculations, Skyrme's interaction would actually 
diverge in perturbation theory because of the zero range. This compels one to use a different interaction to obtain 
correlation beyond the mean field than was used to create it. Extra residual forces are used, such as pairing or 
shell model interactions [ f3l| to allow for more general wavefunctions and obtain more accurate reproduction of the 
physics, or certain approximations are used such as RPA [33]] . which can describe certain observables, particularly 
those of giant resonance states and other forms of collective motion, but not others. It is thought that correlations 
should be particularly important in nuclei at the limits of stability, where the nucleons nearest the Fermi level couple 
strongly to the continuum. 

We revisit the idea that it is possible to parameterize a nuclear interaction in such a way that it is weak enough 
with which to perform perturbation theory, thereby allowing correlated wavefunctions and observables to be calculated 
across the entire range of nuclei. Using the separable ansatz of previous "weak" interactions we have developed a 
density dependent interaction which we hope will provide some insight into the correlation structure of nuclear 
wavefunctions while retaining the quantitative power of contemporary effective interactions used in Hartree-Fock. In 
contrast to previous work, the interaction is designed to be an effective interaction with parameters fitted to the 
properties of finite nuclei within the calculation framework for which it is intended. 

The paper is organized as follows. In Section 2 we briefly review basics of many-body perturbation theory. The 
separable interaction, used in the present work is given and discussed in Section 3. Results of the calculation for 
doubly magic nuclei are summarized in Section 4. Derivation of the HF energy and potential is outlined in Appendix 
A. 



II. SINGLE PARTICLE MODELS AND MANY-BODY PERTURBATION THEORY 



For standard perturbation theory to be successful, the Hamiltonian must be separated into two parts, one which 
is solvable (Hq), and another which is "small" (-Hi). The Hamiltonian for a many-fermion problem in which the 
particles interact via one- and two-body interactions may be written schematically as 

H = X>TO a i a i + \ J2(ij\V\kl)4a] ai a k . (1) 

ij ijkl 

Typically, the one-body part of the Hamiltonian is just the kinetic energy. This is a unsuitable choice for Hq in 
the case of nuclei since the eigenstates of this operator, namely plane wave states, are not close enough to the exact 
eigenstates of the full Hamiltonian and Hi is thus not small. By adding and subtracting a one-body term, U', of ones 
choice, the Hq part of the Hamiltonian may be solvable to give a wavefunction close to the solution of the full problem, 
thereby making the corrections from the Hi part small enough for perturbation theory to succeed. Schematically the 
Hamiltonian is now split up into two terms; 

H = J^(i\U + U'\j) a\a } + i £<ij|U|fcZ)akV* ~ E<^'b'> 4«, ■ ( 2 ) 

ij ijkl ij 
v v ' s , ' 

H Hi 
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One practical choice of U' is a simple, analytically-solvable external potential, typically - in the nuclear case - of 
harmonic oscillator or Woods-Saxon form. In such a case, all of the physics of the nuclear interaction is in the 
residual part, Hi. Unfortunately, for sensible choices of two-body interaction, the residual part of the Hamiltonian 
is non-perturbative in the basis obtained from the solution of the one-body part of the Hamiltonian, and one has to 
perform more exact and exacting interaction shell- model calculations to arrive at a meaningful result. Traditionally, 
this has involved the diagonalisation of large matrices, or more recently, auxiliary field Monte-Carlo calculations in 
the SMMC @. 

Another popular choice for U' is that of the Hartree-Fock mean-field. The HF mean-field potential is usually 
derived from the full Hamiltonian by a variational principle. Viewed in this way, it is the one-body potential whose 
occupied eigenstates form the lowest energy Slater-Determinant many-body wavefunction possible for the full many- 
body Hamiltonian. Using such a one-body potential, some physics of the two-body interaction is included in the 
single-particle problem. With a judicious choice of interaction, one ought to be able to produce very good results, 
since the approximation of the nucleus as a system of non-interacting particles in a mean-field is known to be a good 
one. Using the HF potential for U' has the added attraction that it leads to vanishing first order corrections to the 
energy in perturbation theory. In this case, using Wick's theorem J3f|, one can re- write (0) as 

H = E +^T / e l :a\a t : + ^2{ij\V\kl) :a\a]aia k : (3) 

i<€p ijkl 



H Hi 

where Ho is the HF mean-field Hamiltonian and Hi is the perturbing Hamiltonian. Note that Hi is just the full 
two-body interaction matrix elements with a time-ordered product of creation and annihilation operators. The 
perturbation series for the energy is ordered by the number of matrix elements of the potential. Usually a diagrammatic 
representation is used p6| . One can write down the number of diagrams for any particular order of perturbation theory. 
We have evaluated these for the vacuum amplitude up to order seven. The results are presented in Table |. One sees 
that, for the method of direct evaluation of diagrams by order to be effective, the series must be sufficiently converged 
by fourth, or perhaps fifth, order. 

In this work, the second and third order diagrams for the vacuum amplitude, as shown in Figs. |l| and |^, are 
calculated. Their algebraic form is given here as 

e 2 = 1 V V K°^l">l 2 (4) 

4 £f ^ e a + e b -e r -e s 

ab<t F rs>e F 

E ( PP ) = i y y y {ab\V\rs)(rs\V\tu)(tu\V\ab) 
3 8 ^ ,4" (ta + £b - £ r - e s ){e a + e b ~ e t - e u ) 

ab<e F rs>e F tu>ep v ' x ' 

E (hh) = 1 y y y (ab\V\rs)(cd\V\ab)(rs\V\cd) 

ab<e F cd<ep rs>n F v /x ' 

E ( P h) = y- y- (ab\V\r S )(cr\V\at){st\V\cb) 
3 7^ {e a + e b - e r - e s )(e b + e c - e s ~ e t ) 

abc<e F rst>e F x /v ' 

in which the tildes over the potential indicate the matrix element is antisymmetrized. The state vectors label HF 
single-particle states, whose energies are given by the subscripted e. 

It is important to note that our interaction is not intended to fit scattering data, having, as it does, density 
dependence. On the level of the perturbation theory it is necessary to treat the density functions as just the spatial 
form of the interaction, rather than a representation of a many-body force. This is to be considered a part of the 
present model. To do otherwise would be to surrender the simplifications our weak, separable potential affords. 

In the present work, the Hartree-Fock problem is solved in a basis of spherical harmonic oscillator states. This 
yields, along with the hole states, a large number of particle states, with which to directly evaluate the sums of the 
perturbation series. A sufficient number of states is used so that the particle states are oscillatory over the size of the 
nucleus and that both the HF solution and the perturbation corrections are reasonably converged. 



III. INTERACTION 



We have developed an interaction written in the form of a sum of separable terms, which is to say it is in the form 
Viri, ra) ~ X)s( r i)9( r 2)- The functions g carry no angular momentum (/ = 0), and the force is dubbed a monopole- 
monopole interaction. For future applications, it is in intended to include higher multipole forces, with Z = 1,2,..., 
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within our framework as these will presumably be necessary for calculation of excited states and deformed nuclei. 
Although higher multipole forces will contribute to spherical nuclei from the exchange term in Hartree Fock order 
and via correlations in perturbation theory, they are not included in the present calculation since it seems unwise to 
attempt to fit the parameters of such forces to spherical nuclei alone. 
In coordinate space, the monopole interaction is written as 

V(n,f 2 ) = W a / /-(ri)p A (r2)(l + a a (T+T 2 - + t^t+) + b a t lz t 2z ) 
+ VK r / r p /3 ''(r 1 )p' 3 "(f 2 )(l + a r (T+Tz + t^t+) + b r t lz t 2z ) 

+ k^pin^p^), (8) 

where the function is defined as 



k = 



-l 



<rrp a *(r) 



(9) 



for subscripts £ = a and £ = r. Throughout this work, the three terms in (g|) are referred to, in the order they appear 
in the above expression, as the attractive, repulsive and derivative terms. 
In addition, the spin-orbit force is taken to be 

V s - {r) = c-^l-s, (10) 
r Or 

which is similar to that used in the modified delta interaction J24J . 

The parameters W a , a a , Pa, aa, b a , W r , a r , Pr, cir, b r , k and c are to be fitted to experimental data. 

One notices that the two-body interaction consists of a sum of terms, each of which is separable in form and that 
the expressions for the attractive and repulsive terms in (H) differ only by the values of their parameters 

The ener gy, E pot , due to the interaction (||) in the Hartree- Fock approximation is derived in Appendix |A| QA2| ), 
([A3|), ( A15 ) and is presented here; 



E HF =T + E coul + E pot =T + E coul + {^WtfrNl - \w^M^ + ^W^f^AN^) 2 

£ — a , r 

-\w 6 f 6 [b^M 6 + a $ M^}} + \kNj + cN w (11) 



where T is the kinetic energy, -E C oui is the direct Coulomb energy plus exchange in the Slater approximation. The 
following quantities have been defined; 



N e = J d 3 f/A +1 (r) 

M e = J /d 3 rid 3 f 2 [Mfi^/Hn)/^^ 
AN $ = f d 3 f/«(r) [ Pp (r) - p n (r)} 



= I I d 3 nd 3 f 2 [p p (?u? 2 )pPt(?l)p t(? 2 )p n (?u? 2 )+p n (?l,? 2 )p t^ 
N d = J d 3 rp(r)V 2 p(f) 



N w = / d 6 f-^p w {r) (12) 
J r or 

and the following densities are used: 

P(r) = P P (r) + p n (r) = 'Pi ( ff )'Pi( f ) + Y fii^fiir) (13) 

p(ri,r 2 ) = p p (r u r 2 ) + p n (r u r 2 ) = ^ <Pi {n)<Pi(r 2 ) + Y <Pi (fifyifo) (14) 

P™(f) = H\ + J ) - + V - 3 / 4 ) <Pilf)<Pi<f)- ( 15 ) 
i<t F 
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The variation of the total energy is carried out in Appendix [A| (see |A9] , All, A13, A18) . The resulting local 
Hartree-Fock potential is 



^ — a , r 



- Wtfrft [(1 + bt)G € (x) + a $ G { [ f) (x) 

+ [W i b i P i f i AN i }p^- 1 (x)6p(x)} 
+ 2kN d V 2 p(x) 

which differs for protons (r = p) and neutrons (r = n) through the function 



(,:) 



AN, 



P 

■ n. 



A ^ = { -AN 6 , r 
The other newly-introduced functions in (|l^) are 

G 6 (x) = G\ pp) (x) + G{ nn \x) = J dV [p p (r,x)p^ (r)p p (x,r) + p n (r,x)p^ (r)p n (x,r)} 

G { f\x) = G<f n \x)+G<f p \x) = J d 3 r[p p (f,^(r)p n (x^+ Pn (r,x)p^(^p p (x,r)] 

In addition, the non-local component to the mean field is (see |A14[ ) 

U NLtT (x,x')= W x fJ^{x)^{x') {{\ + b i )p T {x,x l ) + a iP r{S,x)} . 

and there is a state-dependent potential from the spin-orbit interaction of the form 



(16) 



(17) 



U so (x)ip b (x) = c ( w b 



1 dp 1 d Pu 



1 



■ dx 



dx 



Pw{x) (Pb(x) 



(18) 



(19) 



(20) 



where Wb = l/2(jb(j b + 1) — h(h + 1) — 3/4) is the spin-orbit weight factor. 

Note that the one-body spin-orbit term could be taken as either a one-body force, or as a one-body potential 
deriving from a two-body force. Since the latter approach would render the perturbation calculation problematic 
due to the absence of a suitable of form of the two-body force, we choose the former approach. Hence, since the 
force is density-dependent, we have also included the rearrangement contribution to the HF potential. Only the 
non-rearrangement term actually gives rise to the spin-orbit splittings, but the rearrangement terms, coming as they 
do from a variational principle, result in a lowering of the HF energy. Combining the potentials (|l^), (|l^) and ( |20| ) 
gives us the HF equation 



U L , T (x)(pb(x) + / d x' U NL .T(x,x')ip b (x') + U so (x)(f b (x) = Sbfb{x) 



(21) 



In this potential, as well as in the expression for the total energy (|1 1|) , the exchange contribution from the derivative 
term is omitted. While it would, in principle, be desirable to include this term, the calculational complexity involved 
in doing so has forced the omission in the present case. However, for the main attractive and repulsive terms, the 
exchange part is much smaller than the direct in all nuclei, and the direct derivative term gives a rather small 
contribution to the mean-field and the binding energy in comparison to the other direct terms, so it is not considered 
an unwarranted approximation to neglect the effects of this term. 

It might be objected that the form of the interaction is too unrealistic. For instance, since it is separable it cannot 
satisfy Galilean invariance. Furthermore, the unusual form of the isospin operators gives rise to a different effective 
interaction for protons and neutrons. As for the second point, since this is a density-dependent interaction, the effect 
of the difference between neutrons and protons comes from the densities as well as the isospin operators. This being 
the case, it seems reasonable to allow for the "stretching" of the isospin operator as we have done, to allow greater 
freedom in fitting the interaction to data. As for the first issue, the unusual form has been motivated by the desire to 
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keep the force "weak" and is justified by the quality of the results produced. A possible method for restoring Galilean 



invariance within our framework is discussed in Section [V 



The choice of omitting a spin-spin (o ■ a) type force yet having an isospin-isospin type force is motivated by the 
nuclei under study. All the closed-shell nuclei are spin-saturated and would contribute only through the exchange 
term in the HF order. For this separable interaction, the space-exchange terms are rather small and a spin-spin force 
would add little to the results. In addition, even if the effects in closed-shell nuclei are important, it does not seem 
reasonable to fit this term to closed-shell nuclei alone. It remains an open question whether such a force will prove 
necessary or useful in open-shell nuclei. 

It is interesting to compare the leading terms in the HF mean field to that of other models. The first line of equation 
( |l6| ) gives us this as 

where Cj is a combination of constants. The product f^N^ is 

fd 3 r/e +1 (r) 

JK 4 fd 3 rp a i(f) y ' 

If 0£ + 1 = ate then the product /j Ne is constant and the leading mean field terms go like 

U[x) ~ C a pP' (x) + C,./" (x) (24) 

which, for the special case /3 a = 1, are the same as the terms in the Skyrme and Gogny mean-field proportional to 
the parameters to and t^, which give the bulk of the binding energy and saturation properties. In this work, we do 
not strictly keep (3^ + 1 = thus allowing for some A-dependence of the coefficients in the mean- field potential. It 
has been found, however, that one can not go too far away from the equality and still obtain reasonable results. 

That one can get similar results in a mean-field calculation from two very different interactions is reflected in the 
different constitution of the resulting perturbative part of the Hamiltonian, Hi. 



IV. DOUBLY (SEMI-) MAGIC NUCLEI 

In order to find the best set of parameters for the interaction (|J), calculations have been made of 14 doubly closed- 
shell nuclei across the periodic table. They are 16 0, 34 Si, 40 < 48 Ca, 48,56,68,78 Ni; 90^ ioo,iu,132 Si1j i46 Gd and 208 pb 
The nuclei represent a selection of doubly-closed (sub-) shell nuclei both close to and far from stability. There is limited 
experimental information about 48 Ni []42| and 100 Sn ||3|. 78 Ni has yet to be discovered. 

The ability to reproduce the properties of such exotic nuclei will be important for applications of our technique and 
discrepancies will help direct refinements. 

A Hartree-Fock code assuming spherical symmetry and representing wavefunctions in a basis of spherical harmonic 
oscillator states was used to calculate uncorrelated wavefunctions. Perturbation corrections to the binding energy 
were directly evaluated using the results of the HF calculation. The results presented here were obtained in a basis of 
12 expansion coefficients per single-particle wavefunction and iterated until the HF energy had converged to within 
IkeV. The parameters of the force were fitted to binding energies to second order and charge radii, charge density 
distributions, single-particle energies and spin-orbit splittings to HF order of the nuclei listed above where experimental 
data were available, and are presented in Tabic ||. 

The results of the the calculated energies, in HF order and in each order of perturbation theory are presented in 



Table |Hj. The differences between the HF energy E ca \ = (11|)-Bhf an d the experimental ground state energy i? cxp and 
between HF plus second order perturbation correction (^ ) E ca i = Ehf + Ei and experiment are shown in Figure |^. 
The experimental energies are taken from the mass table of Audi and Wapstra Gq] with two exceptions. An estimate 
of the mass of the recently-discovered nuclei 48 Ni (5l) and the measured mass of™Sn j|4|. The energy for 78 Ni was 
taken from |45| ] in which extrapolated values are given, which are thought to be in error by less than 0.2%. 

One sees from Figure || that most of the nuclei fit the binding energy to within ~2%. The most obvious exception 
is 16 which is quite under-bound. This may be due to the omission of a center-of-mass correction which would 
undoubtedly go a long way to close the discrepancy in the energy |4(J. It was not calculated in this case since a 
rigorous microscopic correction would destroy the mean-field which provides the essential basis for the perturbation 
calculation. A phenomenological correction could have been calculated, but perhaps the most suitable method in our 
framework would be to include with our multipole forces an isoscalar dipole force which could be fitted to restore the 
translational invariance of the many-body Hamiltonian, and evaluated exactly in perturbation theory. 
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A general trend can be seen in which lighter nuclei are somewhat over-bound and the heaviest are under-bound. It 
is the exceptions which conspire to stop the fitting algorithm from doing better, but the somewhat systematic nature 
of this discrepancy suggests that a better mass or isospin dependence may improve matters. It is unclear as yet the 
extent to which multipole correlations or a spin-spin force would improve the fit to spherical nuclei. That question 
awaits the study of deformed nuclei and excited states. 

In Table |y] a comparison is made of the quality of the fit to the binding energy to properties of the same nuclei 
calculated with a selection of Skyrme parameterizations. The parameterizations used are SIII jl?]], SkP [jl8|, SLy4 
p9| and SkI4 f30| . In this comparison, it is seen that the energies from the different Skyrme parameterizations are 
of a similar quality, all reproducing the binding energies of closed-shell nuclei very well, with only a few binding 
energies being reproduced no better than 1% - including 48 Ni whose experimental value is in any case not well known. 
It is clear that the results from the separable force are somewhat worse. Particularly problematic is 16 O, whose 
large under-binding was mentioned above, and also 48 Ni which is, as with the Skyrme parameterizations, over-bound, 
although more so with the separable interaction. 

Results for one-body properties are also presented. Comparison of the charge density results to experiment |5q] 
and to the selection of Skyrme interactions is made in Table [vj]. One-body observables are generally reproduced 
better in the HF calculation alone than the binding energies. The comparison of the calculated charge radii with 
experiment is generally more favourable than the energy data. The radius of oxygen is too large by about 5% which 
is consistent with its under-binding. It can be seen that the agreement with experiment is of the same level as the 
Skyrme interactions. Perturbative corrections to the one-body observables, such as the densities, and hence radii, will 
be calculated in future work. 

Figures |^-|| show the electron scattering form-factors of a selection of the nuclei compared to experiment |55| . The 
proton density was corrected for the finite proton size by folding with a Gaussian to give the charge density, from 
which the radii and form-factors were calculated. The form-factors agree with experiment quite well, which is expected 
given the generally correct radii. 

Some typical single-particle energies for light and heavy nuclei are shown in Figures and |l] for 40 Ca and 208 Pb 
respectively. The single-particle energies of a density-dependent Hartree-Fock calculation do not directly correspond 
to an experimental observable, so caution should be used in comparing values. It can be seen that the level spacings 
and shell closures are better reproduced in 208 Pb, which is true of heavy nuclei in general. The comparatively poorer 
results in light nuclei seems to be common to mean- field approaches [ pl| . In the case of 40 Ca, the gap at the fermi 
level may be widened by the inclusion of multipole forces which will link the occupied levels with the states, 
which is otherwise "inert" . 

Table shows spin-orbit splittings for some cases where the experimental values are known. The "experimental" 
data presented represents that used in previous work for fitting effective interactions to data ^||56|,^7j . Clearly the 
splittings are all systematically small. This could be remedied by an increase in the spin-orbit coefficient, c. In a 
previous work (5£| , a value 10% higher than ours was used for the same spin-orbit interaction, and hence the spin-orbit 
splittings were more realistic. The lower value used in our work is the result of a compromise between the reproduction 
of the spin-orbit splittings and the total binding energies. This is a further indication that a more suitable spin-orbit 
potential needs to be sought. 

The perturbation correction to the energy are seen to be rather small in all nuclei considered. This is consistent with 
our goal that the mean field solution should be close to the exact solution of the MBSE. The size of the second order 
correlation is roughly constant across the periodic table. It is characterized by a dimensionless strength parameter, 
k, defined as 

k= 1_ E E (ab\V\rs)(rs\VW ^ 

It is related to the "wound integral" |5j] and is proportional to the number of 2p2h states excited due to the second 
order perturbation. 

Figure shows the correlation structure from the second order correction in the N = Z nucleus 40 Ca, and the nuclei 
48 Ca and Pb. The contribution to the second-order energy is defined as a function of one of the particle states r; 



i \{ab\V\rs) 

ab<£F S>£F 

The plot shows the contribution to the total second order energy correction as a function of the single particle energy 
e r , in 5 Mev wide bins. In all three cases particles are dominantly excited to low-lying states above the Fermi level. 
This results in a ground state with occupation probabilities similar to those which result from pairing forces. It is 
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also a further indication that perturbation theory makes sense for our interaction since it does not predict excitation 
of particles in the ground state to extremely high energies. 

Since the calculations were made only using a monopole force, the correlation structure is not expected to be 
complete. Only corrections involving simultaneous I = scattering of two particles is included. An indication of this 
is seen in the difference between the results for N = Z and N ^ Z nuclei. The second order correction in 48 Ca is much 
larger than that in 40 Ca, due to the possibility of an f^/ 2 neutron exciting to the proton state while another proton 
excites to a neutron state. This extra excitation is the labelled peak in Fig. |[ As well as having large wavefunction 
overlaps, the energy denominator in this case is much smaller than in any other possible monopole excitation which 
must excite any particle across major shells to keep all angular and isospin quantum numbers the same. When general 
excitations are permitted by higher multipole forces I = 1,2..., this difference between correlations in N — Z and 
N ^ Z nuclei will be smoothed out. For this reason, too, the correlation energies should not be considered too 
quantitatively at this stage, but rather as an indication of the perturbative properties of the interaction. 

Comparing the form of the interaction to that of Skyrme suggests other possible sources of improvement to the 
model. One such may come from a better parameterization of the spin-orbit interaction. A two-body form which 
fits the philosophy of the separable effective interaction has not been found, but may be necessary to give the correct 
contribution to the binding energy. In any case, the simple form used in the present work which depends on a radial 
derivative will need modification if it is to be applied to deformed nuclei. It may also prove fruitful to explore a more 
general term dependent upon the derivatives of the density than the single term with parameter fc, such as is found 
in the Skyrme interaction with two terms proportional to ti and t 2 , which often carry further exchange parameters 
x\ and X2- 

V. CONCLUSION 

We have presented a new effective nuclear interaction which is designed for use in calculations which go beyond 
the mean-field. The technique of using perturbation theory to build correlations on top of the Hartree-Fock result is 
applicable to our interaction and results in small corrections to the single-particle behavior. A monopole- monopole 
force alone gives reasonable results for the ground state properties of spherical doubly-magic nuclei. It is expected 
that the addition of multipole forces will improve these results, particularly through the completion of the correlation 
structure. Such multipole forces will also presumably be important in giving the correct shapes of deformed nuclei, 
which are the subject of a forthcoming study and in the correct description of excited states. 
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APPENDIX A: HF ENERGY AND POTENTIAL 



The interaction is given in equation (M). Its expectation value is the contribution it makes to the total energy and 



is 



E V ot = \ £ (ij\V{\ij)-\ji)}. 



(Al) 



ij<e F 



If we consider just the attractive term - that is, the term whose parameters have the subscript "a" - then we will 
obtain the contribution from the repulsive term by simply substituting the subscript "a" for "r" . 

E a = \wj a ]T (ii|^(n)/«(r 2 ) - \ji)} 

ij<t F 

+ \w a f a a a fo> & (ri)/°(r 2 )(T+T 2 - +rf r+){\ij) \ji)} 

ij<£F 

+ \w a f a b a (ij\p 0a in)p^(r2)ti z t 2z {\ij)-\ji)}. 



(A2) 



ij<a F 



Taking the first line, the matrix element is represented in space (and spin and isospin) coordinates; 

\w a f a J dV i / d 3 r 2 rAri)<f>*(r2)p Pa (n)p Pa (r2)Mri)Mr2) 



ij<e F 



\Waf a J dV l / dV 2 ^(ri)^(f 2 )p A (fi)p^(f 2 )^(fi)0 i (f 2 ) 



ij<a F 



\Waf a ( dM A +1 (rl) / dV 2 p^ + \r 2 ) 



~W a f a I dM / d a r 2 Pv {f u i t 2)p^{fi)pP*{f2)Pp(f2,f , i) 



W a f a I dM / d 3 r 2 Pn (r 1 ,r 2 )p^(r 1 )p^(f 2 )p n (r 2 ,r 1 ) 



\Waf a N 2 a - ~W a f a M a 



(A3) 



where quantities defined in section III are used. Note that integrals include sums over spinors and isospinors where 
appropriate, and the coordinates include spin and isospin coordinates where appropriate. Where densities are used, 
the summing over isospin states has already been done and where densities do not carry isospin labels, the isoscalar 
density is assumed. S ee e quations (13) and (f[4|) for definitions. 

The second line of (A2) contains "isospin-flipping" operators whose action is to turn an isospin state where particle 
one is a proton and particle two a neutron, \pn), into \np) and vice- versa. The direct contribution, in which the labels 
in the bra and the ket are in the sa me order is zero since all proton states are orthogonal to all neutron states. The 
exchange term is similar to that in (A3) but with different isospin combination of the density matrices; 



~W a f a a a I d a ri / d 3 r 2 p p (f u r 2 )/A (ri)p Pa (f 2 )p„(f 2 , fi) 



--^W a f a a a j dVi / dV 2 Pn(ri,f 2 )p f>a (fi)p^(f 2 )pp(f 2 ,f 1 ) 
= -\w a f a a a M^ 



(A4) 



The third line of A2 contains isospin-projection operators which have a value +1 when i and j are like particles and 
-1 when they are unlike. In the direct term this gives an energy of 

Wafaba [ d 3 fi ft,(fi)/° (n) / d 3 r 2 ft, fa)/" (r 2 ) 







+ ^W a f a b a J dft^ftV^ft) J dV 2Pn (r 2 )/»(r 2 ) 

-\Wafaba J d 3 f=i P p ftl)/Aftl) ^ dV 2 p„ ft )p A ft ) 

1 



W Q / a 6 a / d 3 fi/> n ft)p^ft) I d A r 2 p p {r 2 )p^{r 2 ) 



^W a f a b a 



d 3 rp p (r)p^(r) 



^W a f a b a 



d 3 rp n (r)p f) "(r) 



1 



W a f a b a / d 3 ^ Pn^l)/^' ft) / dft ft, (5)/° ft) 

n 2 

d 3 f( Pp (r)-p„(r))p^(r) 



= ^W a f a b a 



1 



|W«/ &«(AJV a ) a . 



(A5) 



The exchange term gives a contribution only when i and j have the same isospin quantum number, which is just like 
the case for having no isospin operator there, so the contribution is like that in (A3) 



\w a f a b a M a . 



(A6) 



The HF mean-field is obtained by varying the total energy with respect to the single particle states. This gives, for 
the case of the attractive term, without the explicit isospin dependence: 



u b (x) \2 2 



l W Sfa N 2 + W f N SNa 

l w Sf a ^ 1 6M a 
2 5<t>l (x) 



The variation of the function /„ is given by 



Sfa 


6 




S4>t(x) 


S(f>l(x) 


:/ 




-fa J d»f 




~fa a a J 





p aa {r)d 3 r 
tyj (2) 

fla a p^- x {x)()> b {S) 



(A7) 



, Sp(f) 



(A8) 



so that the contributions of the two terms in ( A7) involving the variation of f a give a contribution to the HF mean field 
of 



-W a {a a /2)ft{Nt-M a )p<*«-\x) 
The function N a is similar in form to f a and the functional variation proceeds in a similar manner; 

SN a S f 3 _ g 

d 3 r p Pa+ \r) 



(A9) 



= O8 o -l)/«(f)0 6 (5) 



and the contribution from the second term in (A7) to the mean-field is 

W a f a N a {p a -l)pP«{x). 
Finally, the functional variation of the exchange matrix element M a is 



(A10) 
(All) 
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SMg 6 



J2 dMdV^^n)^^)^^)^^)^^)^^) 



lJ<tF 



2 E / /d 3 rldVl^^0^f 2 )p^(f 1 )p^(^)0i(n)^(^) 



2 £ //d^d^anj^ft)^^^),,!^,,.^) 



lJ<t F 



2 I d 3 r 2 0*(f 2 )p^(x)p^(r 2 )^(f)0 fc (f 2 ) 

ij<e F 



(A12) 



where use is made of the symmetry of the integral to combine the four terms into two. The last term in (A12) gives 
rise to a local term in the mean-field of 



i<£F 



-W a f a f3 a G a (x)p^'\x) 



where G a (x) has bee n defi ned as in (|18[). 

The other term in (A12) gives rise to a truly non-local Fock term in the mean-field: 

U(x,x')Mx) = -W a fa 53 J d 3 f^*(rV^(r)<Mr) 



i<£F 



(A13) 



(A14) 



This completes the non-isospin-dependent part of the attractive force, and so also the repulsive by change of 
subscript. The isospin-dependent terms are obtained in an analogous way, except that when the variation applies to 
the density of a single nucleon species, so the contribution to the mean-field applies only to that species. 

For the final term in equation (^), the so-called derivative term, only the direct part of the energy is at present 
considered. It is 



£dcriv. = \k J2 fe'|Vi 2 p(^)V 2 2 (r 2 )|ij) 



k / dMp(riKp(ri) / d d f 2 p(r 2 )\7 2 p(f 2 ) = -kNl 



(A15) 



The functional variation proceeds as 



SE dc 



kN rt 



d 3 fp(f)W 2 p(f) 
Sp(r) 



kNd I dV imk } yW) + kNd I dV>(r1 {sm y2p(f) 



The first term gives a contribution to the mean field of 

kN d V 2 p{x). 



(A16) 



(A17) 



By integrating the second term by parts twice, one in fact gets exactly the same contribution to the mean-field again, 
so that the total contribution to the mean-field from the direct term of the derivative interaction is 



2kN d V 2 p(f) 



(A18) 



The exchange part of this term is not calculated. 
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TABLE I. Number of Hugenholz diagrams by order of perturbation theory 



Order 


2 


3 


4 5 


6 


7 




Labelled Diagrams 
Cumulative Diagrams 


1 
1 


3 
4 


39 840 
43 883 


27,300 
28,183 


1,232,280 
1,260,463 








TABLE II. 


Monopole force parameters 








W„ 

-1543.8 MeV fm 3 




- * a 

2.0 


B„ 
1.0 


n - 

-0.4295 




b„ 

-0.444825 


W r 

1778.0 MeV fm 3 - 8265 




a r 
2.2165 


Pr 

1.246 


a r 
-1.4788 




b r 

-0.314625 


c 

160.0 Mev fm 5 








fe 

16.0 Mev fm 10 







TABLE III. Hartree-Fock energy ( |ll[ ) and corrections from perturbation theory (y)-(0) compared with experimental value 
from All energies are in MeV 



Nucleus 


Ehf 






^(3) pp 




Eh f+2+3 


K 


Expt. 


ib 


-109.32 


-3.31 


-0.1365 


-0.3624 


+0.921 


-112.21 


0.063 


-127.68 


34 Si 


-280.88 


-7.37 


-0.0384 


-0.4830 


+ 1.223 


-287.55 


0.232 


-283.43 


40 Ca 


-334.53 


-2.51 


-0.0323 


-0.1114 


+0.233 


-336.95 


0.052 


-342.00 


48 Ca 


-417.01 


-5.97 


-0.0189 


-0.2725 


+0.273 


-422.70 


0.202 


-416.16 


48 Ni 


-360.69 


-6.57 


-0.0130 


-0.2058 


+0.427 


-367.05 


0.234 


-348.33 


56 Ni 


-481.25 


-2.31 


-0.0210 


-0.0643 


+0.123 


-483.52 


0.046 


-483.99 


68 Ni 


-593.33 


-6.00 


-0.0109 


-0.2091 


+0.484 


-598.85 


0.221 


-590.43 


78 Ni 


-651.90 


-8.34 


-0.0053 


-0.1458 


+0.477 


-659.92 


0.342 


-641.38 


90 Zr 


-782.70 


-3.91 


-0.0070 


-0.1257 


+0.103 


-786.51 


0.149 


-783.89 


100 Sn 


-825.65 


-1.71 


-0.0060 


-0.0220 


+0.048 


-827.35 


0.039 


-826.81 


114 Sn 


-963.20 


-4.04 


-0.0046 


-0.1093 


+0.226 


-967.12 


0.162 


-971.57 


132 Sn 


-1097.65 


-6.17 


-0.0023 


-0.0864 


+0.209 


-1103.70 


0.287 


-1102.92 


146 Gd 


-1190.32 


-3.42 


-0.0026 


-0.0699 


+0.142 


-1193.66 


0.146 


-1204.44 


208 pb 


-1599.04 


-4.51 


-0.0013 


-0.0664 


+0.108 


-1603.51 


0.233 


-1636.45 



TABLE IV. Percentage error in binding energy. Negative values are underbound. Separable force is HF+perturbations, 
Skyrme calculations are HF+pairing. 



Nucleus 


Sep. 


SIII 


SkP 


SLy4 


SkI4 


16 


-12.10 


0.36 


-0.12 


0.19 


0.57 


34 Si 


1.45 


0.43 


0.92 


1.08 


1.06 


40 Ca 


-1.49 


-0.18 


0.21 


0.50 


0.51 


48 Ca 


1.61 


0.40 


-0.04 


-0.63 


0.31 


48 Ni 


5.37 


1.52 


1.22 


0.74 


1.58 


56 Ni 


-0.10 


-0.22 


-1.11 


0.09 


-0.29 


68 Ni 


1.43 


-0.25 


0.09 


0.81 


0.26 


78 Ni 


2.89 


0.58 


-0.05 


-0.29 


0.21 


90 Zr 


0.33 


-0.14 


-0.14 


-0.11 


0.13 


100 Sn 


0.07 


0.14 


-0.51 


0.75 


0.28 


114 Sn 


-0.46 


-0.71 


-0.48 


0.09 


-0.56 


132 Sn 


0.07 


0.10 


-0.31 


0.14 


-0.12 


146 Gd 


-0.90 


-0.33 


-0.41 


-0.20 


-0.24 


208 pb 


-2.01 


-0.17 


-0.27 


0.21 


-0.24 
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TABLE V. Spin-Orbit splittings in HF calculation. For source of experimental values, see text. 



levels 



O, Op (n) 



40 Ca, Od (p) 
40 Ca. Od (n) 
208 Pb, 2p (n) 



splitting (HF) 

u~ 

4.3 
5.3 
5.3 
0.67 



splitting (exp) 
63 
6.1 
7.2 
6.3 
0.89 



TABLE VI. Comparison of charge radii between experiment, the separable interaction and a selection of Skyrme interactions. 
The model-dependent experimental values are from Eq{ 



Nucleus 


exp. 


Sep. 


SIII 


SkP 


SLy4 


SkI4 


16 o 


2.69 


2.85 


2.71 


2.80 


2.76 


2.72 


34 Si 




3.19 


3.23 


3.25 


3.23 


3.21 


40 Ca 


3.48 


3.54 


3.48 


3.52 


3.49 


3.45 


48 Ca 


3.48 


3.47 


3.52 


3.53 


3.51 


3.45 


48 Ni 




3.88 


3.77 


3.82 


3.79 


3.80 


56 Ni 


3.78 


3.84 


3.80 


3.80 


3.78 


3.74 


68 Ni 




3.89 


3.94 


3.93 


3.91 


3.81 


78 Ni 




3.87 


4.02 


3.99 


3.98 


3.97 


90 Zr 


4.27 


4.29 


4.31 


4.30 


4.28 


4.23 


100 Sn 




4.60 


4.53 


4.52 


4.50 


4.45 


114 Sn 


4.60 


4.65 


4.66 


4.62 


4.62 


4.59 


132 Sn 




4.66 


4.78 


4.74 


4.73 


4.70 


146 Gd 


4.96 


5.02 


5.03 


5.00 


4.99 


4.94 


208 pb 


5.50 


5.50 


5.57 


5.52 


5.51 


5.48 
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FIG. 1. Labelled second-order Hugenholz diagram 




FIG. 2. Labelled third-order Hugenholz diagrams 
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200 



FIG. 3. Deviation of Hartree-Fock Energy from experiment. Negative errors denote underbinding. Note that le O and the 
second-order result for 48 Ni are beyond the scale. 
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q [feT 1 ] 

FIG. 5. Charge Form factor in 16 O 
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q [fnf '] 

FIG. 8. Charge Form factor in 90 Zr 




q [mf '] 

FIG. 9. Charge Form factor in 208 Pb 
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FIG. 10. Single-particle energies in Ca compared to experiment |51 




FIG. 11. Single-particle energies in 208 Pb compared to experiment 
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